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We investigate the influence of different interaction strengths and dimerizations on the magne- 
tization transport in antiferromagnetic spin 1/2 XXZ-chains. We focus on the real-time evolution 
of the inhomogeneous initial state | | . . . || . . . |) in using the adaptive time-dependent density- 
matrix renormalization group (adaptive t-DMRG). Time-scales accessible to us are of the order of 
100 units of time measured in h/J for almost negligible error in the observables. We find ballis- 
tic magnetization transport for small aS^aS^ -interaction and arbitrary dimerization, but almost no 
transport for stronger iS'^S'^-interaction, with a sharp crossover at = 1. Additionally, we perform 
a detailed analysis of the error made by the adaptive time-dependent DMRG using the fact that the 
evolution in the XX-u\ode\ is known exactly. We find that the error at small times is dominated by 
the error made by the Trotter decomposition whereas for longer times the DMRG truncation error 
becomes the most important, with a very sharp crossover at some "runaway" time. Overall, errors 
are extremely small before the "runaway" time. 



I. INTRODUCTION 

The transport properties of spin chains have attracted 
much attention recently, not only due to the possible ap- 
plications to information storage, spintronics, and quan- 
tum information processing, but also because they allow 
to study general aspects of nonequilibrium dynamics in a 
comparably simple system. Nonequilibrium phenomena 
are a vast and despite all progress still poorly under- 
stood field of statistical physics. It is therefore useful to 
have a simple model at hand that allows to study general 
questions rather explicitly. In order to study nonequilib- 
rium phenomena, a real-time description is particularly 
intuitive and useful. In this paper, we study the time- 
evolution of a spin-1/2 chain by solving the full many- 
body Schrodinger equation. 

Recently, new developments in the area of non- 
equilibrium physics were stimulated by the experimental 
progress in the field of ultracold atoms. These systems 
have the advantage that their parameters can be tuned in 
time with high accuracy and on very short time-scales. It 
was proposed that spin-1/2 chains can be realized in these 
systems as well^'^'^'^, namely as a mixture of atoms of two 
species, say A and B. If these atoms are studied in an 
optical lattice with an average filling of one atom per site 
and with a very strong repulsive interaction between the 
atoms, such that multiple occupancy is suppressed, the 
system can be mapped onto an effective spin-1/2 model. 
In this effective model the state with atom A occupying 
a given lattice site corresponds to, say, |, and likewise B 
to j. 

In this paper, we study the time evolution of an initial 
state I t . . . Ti • • • i) (or \A"'AB"'B)), i.e. with all 
spins on the left half pointing up along the z axis, and 
all spins on the right half pointing down, under the ef- 
fect of a nearest-neighbour spin interaction (see Eq. QJ. 
This system can also be interpreted as an oversimplified 
picture for spin transport between two coupled reservoirs 



of completely polarized spins of opposite direction in the 
two reservoirs. We are mainly interested in the follow- 
ing questions: Does the state evolve into a simple long 
time limit? If so, how is this limit reached? On what 
properties does the long time behaviour depend? 

Analytical results for this problem are essentially re- 
stricted to the XX-chain which is amenable to an exact 
solution^. Here, a scaling relation for the long-time limit 
was found. However, it is presently not known whether 
this relation is general, or whether it relies on special 
properties of the XX model. If a long-time limit exists 
for other models as well, the question arises which of its 
characteristics are universal, and which depend on cer- 
tain system properties. 

Directly solving the time-dependent Schrodinger equa- 
tion for interacting many-body systems is highly nontriv- 
ial. A recently developed numerical method, the adaptive 
time-dependent DMRG^'^'^ (adaptive t-DMRG), enables 
us to perform this task. The two main conditions for this 
method to be applicable, namely that the system must be 
one-dimensional and have nearest-neighbour interactions 
only, are met for the present model. Efforts to generalize 
the DMRG method to time-depending problems relaxing 
these constraints are under wa}iS. 

As so far no detailed error analysis of this new method 
has been performed, an important aspect of the present 
work is that besides their own physical interest, spin-1/2 
chains provide an excellent benchmark for the adaptive 
time-dependent DMRG, because of the nontrivial exact 
solution for the XX model, against which the method 
can be compared. This allows us to analyze the accuracy 
of the adaptive time-dependent DMRG very explicitly, 
namely to address the questions what kinds of errors can 
occur in principle, which ones of these dominate in prac- 
tice, and how they can be minimized. We find that the 
time-scales accessible to us are about lOOfi/J with an 
neglegible error in the observables at very moderate nu- 
merical cost. 
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FIG. 1: Quantum phase diagram of the Heisenberg model, 
Eq. (HI. See |ia,llj for details. 



The outline of our paper is as follows: In section Hll we 
introduce the model and its characteristics. In section lTTTI 
we summarize the method and in section Hvl a detailed 
error analysis is performed. These two sections may be 
skipped by readers mainly interested in the physics and 
not in the details of the method. In section Owe present 
our results for the long time limit of the time evolution 
for different interaction and dimerization strength. 



II. 



MODEL AND INITIAL STATE 



In this paper we analyze the dynamics of the inho- 
mogeneous initial state |ini) = | | ... || ... |) on the 
one-dimensional spin- 1/2 chains with interactions given 
by the Heisenberg model 

n n 

Here, Sn is the spin operator on site n, and J^, Jz are in- 
teraction constants. We consider dimerized models where 
= const and Jn = (1 + (—1)^^), S being the dimer- 
ization coefficient. For (5 > 0, the "strong bond" with 
Jn = 1 + is chosen to be at the center, where the spin 
flip of the initial state is located. 

We have chosen our energy unit such that = 1 for 
the homogeneous case S = 0. We also set h = 1, defining 
time to be 1/energy with the energy unit chosen as just 
mentioned. 

The quantum phase diagram of this model at zero tem- 
perature is well known (seeiSiii) and sketched in Fig. ^ 
For the homogeneous case, ^ = 0, the ground state 
has ferromagnetic (FM) / antiferromagnetic (AFM) or- 
der with a gap in the excitation spectrum for < — 1 
and > 1, respectively. The gap closes if \Jz\ ap- 
proaches 1 from above, and the model becomes critical 
for — 1 < < 1, i.e. gapless in the thermodynamic limit, 
with correlation functions showing a power-law decay. 
The model at the point Jz = S = is known as the XX 
model. It has the special property that the spin-current 
operator J = ^^jn is conserved, i.e. [J^H] = 0. Here 
jn = Jnim {S^S~_^i) is the current operator on the bond 
between site n and n + 1. For finite dimerization, ^ 7^ 0, 
the spectrum is again gapped for all values of J^. 



Often it is useful to map the Heisenberg model onto a 
model of spinless fermions: 



-\-Jz{cl^C^ ^) {cl^-^lC^_^l ^) 



(2) 



In this picture, the first two terms in Eq. describe 
nearest-neighbour hopping, whereas the third term (the 
one proportional to Jz) describes a density-density in- 
teraction between nearest neighbours. In particular, the 
case Jz = describes free fermions on a lattice, and can 
be solved exactly 

The time-evolution under the influence of a time- 
independent Hamiltonian H as in Eq. is given by: 



IV^(t)) = U{t)\mi) with U{t) = exp{-iHt). (3) 

In most of the phases shown in Fig. ^ the state 
|ini) = |t...||...|) contains many high-energy excita- 
tions and is thus far from equilibrium. In the following, 
we briefly discuss these phases separately. 

- Deep in the ferromagnetic phase, Jz < —1, |ini) corre- 
sponds to a state with one domain wall between the two 
degenerate ground states. For ^ —00 it is identical 
to the ground state (with boundary conditions given by 
I t) and I I) and 5*°* = 0), and therefore stationary. For 
finite J^, it is no longer identical to the ground state, but 
still close to it^^. 

- In the antiferromagnetic phase, > 1, the state |ini) 
is highly excited. One could view it as a state with al- 
most the maximum number of domain walls of staggered 
magnetization. 

In this context, it is interesting to note that the sign 
of Jz does not matter for the time evolution of physical 
quantities, as long as the initial state is described by a 
purely real wave function (which is the case for our choice 
of |ini)), since the sign change in Jz can be compensated 
by a gauge transformation that inverts the sign of the 
hopping terms S^S^^ S^S^ in Eq. (QJ, plus a complex 
conjugation of Eq. 0. In particular, the time-evolution 
of the low-energy one domain- wall state in the FM is the 
same as the evolution of the high-energy many domain- 
walls state in the AFM. We therefore restrict ourselves 
to the case Jz > 0. 

- In the critical phase ^ = and \Jz\ < 1, the ground 
state is a state with power-law correlations in the xy- 
plane. Here, the state |ini) is not close to any particu- 
lar eigenstate of the system, but contains many excited 
states throughout the energy spectrum, depending on the 
value of Jz'. The energy expectation value of |ini) is low 
as ^ — 1 and high as ^ 1. 

The time evolution delocalizes the domain wall over the 
entire chain. For = 0, the time-evolution of the system 
can be solved exactly. For example, the magnetization 
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FIG. 2: Illustration of the DMRG bases with active site n 
and n + 1, respectively. 



profile for the initial state |ini) reads^: 



s^t) = imis'nim) = -u^ E -^iW' (4) 

j=l-n 

where Jj is the Bessel function of the first kind, n = 
. . . , — 3, — 2, — 1, 0, 1, 2, 3, . . . labels chain sites with the 
convention that the first site in the right half of the chain 
has label n = 1. 

- In the dimerized phase, 7^ 0, the mentioned charac- 
teristics remain unchanged. However, here the delocal- 
ization becomes confined to pairs of neighbouring sites in 
the limit S ^ 1. 

We finally note that the total energy and magnetiza- 
tion of the system are conserved at all times, such that 
even for long times the state cannot relax to the ground 
state. 



III. OUTLINE OF THE ADAPTIVE 
TIME-DEPENDENT DMRG FOR SPIN CHAINS 

In order to determine the time-evolution of Eq. ((H)), we 
use the adaptive t-DMRG method^-^, which has been 
introduced as an extension of standard DMRG using 
the TEBD algorithm of Vidal&. It allows to evaluate 
the time-evolution for one-dimensional quantum chains 
with nearest-neighbour (possibly time-dependent) inter- 
actions. In this paper, we consider the case of a time- 
independent Hamiltonian where the dynamics is intro- 
duced by a nonequilibrium initial state at t = 0. To set 
the stage for the error analysis, we briefly review adaptive 
t-DMRG, assuming the reader to be familiar with stan- 
dard static zero temperature DMRG (see, e.g. 14 15]). 

In the standard finite-system DMRG algorithm, a 
quantum-mechanical state on a one-dimensional chain 
with L sites is represented in a particular tensor prod- 
uct basis, namely as 



OiarjS 



aT(3\(y^)l..n-l\(y)n\Vn+l\p)n+2..L 



(5) 



as illustrated in the upper part of Fig. [21 Here, |cr)n, 
|r)n+i are complete bases on sites n, n + 1; |a)i..n-i and 
|/5)n+2..L are states on the subchains with sites 1, n—1 
and n + 2,...,L, respectively. The states |Q^)i..n-i and 
|/^)n+2..L form truncated bases, i.e. they do not span 



the full Hilbert space on their respective subchains, but 
only a subspace of dimension m, chosen to allow an op- 
timal approximation of the true physical state. In the 
representation of Eq. (O, we call site n the "active site". 
The algorithm now consists of moving ("sweeping") the 
position of the active site several times from the left to 
the right end of the chain and back, and constructing 
optimized truncated bases for the subchains. 

A DMRG step during such a sweep, say, to the right 
now consists of a basis transformation from the old (trun- 
cated) basis \(y)i..n-iW)n\^)n+i\P)n+2..L with active site 
n to anew one \a')i,,n\cr')n+i\r')n+2W)n+2...L with active 
site n + 1 as shown in Fig. [21 The states |a')i..n repre- 
senting the sites 1, n are linear combinations of the old 
basis vectors \(y)i..n-iW)n' Not all linear combinations 
are kept because of the DMRG truncation that limits the 
number of states to m states. For this reason, the 

state 1?/^) can in general be represented in the new basis 
only up to some truncation error. The DMRG trunca- 
tion algorithm (described in |3^3) provides a unique 
optimal choice for the states \a') that minimizes this er- 
ror (which is then typically as low as 10~^^ or so) and 
thus allows for the optimal representation of particular 
"target" states. The basis vectors |/3')n+3..L are taken 
from stored values from the previous sweep to the left. 
- A sweep to the left (i.e. from active site n to n — 1) 
works in the same way, with the role of \a') and 
interchanged. 

In standard DMRG, a mere transformation of the state 
from one basis to the other - known as White's state 
prediction^^- is possible and accurate up to the (small) 
truncation error. However, in order to optimize the basis 
states iteratively for representing the target state (s) |?/;), 
new information must be provided about |?/;), i.e. it must 
be newly constructed using some unique criterion (typi- 
cally as the ground state of some Hamiltonian) . Without 
such a criterion to "sweep against", the accuracy can- 
not increase during sweeps, and the procedure would be 
pointless. Merely transforming |?/^) in this way is there- 
fore of no use in standard DMRG, and is in fact never 
performed alone. It is, however, the basis of the adaptive 
t-DMRG. 

The adaptive t-DMRG algorithm relies on the Trot- 
ter decomposition of the time evolution operator U{t) 
of Eq. which is defined as follows: Using the rela- 
tion U{t) = U{dt = t/M)^, the time evolution operator 
is decomposed into M time steps, where M is a large 
number such that the time interval dt = t/M is small 
compared to the physical time scales of the model. Since 
the Hamilton operator of Eq. (^) can be decomposed into 
a sum of local terms hn that live only on sites n and n + 1, 
U (dt) can then be approximated by an n-th order Trotter 
decompositionii, e.g. to second order: 



U{dt) = l[uAl[Un{dt) l[Un{^] 



■0{dt^). (6) 



odd 



The Un{dt) are the infinitesimal time-evolution opera- 
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tors exp{—ihndt) on the bonds n (even or odd). The 
ordering within the even and odd products does not mat- 
ter, because "even" and "odd" operators commute among 
themselves. 

Eq. ahows to decompose the time-evolution opera- 
tor U (t) into many local operators Un that live on sites n 
and n + 1. The adaptive time-dependent DMRG now al- 
lows to apply the operators Un successively to some state 
^. Each operator Un is applied exactly during sweeps in 
the DMRG step with n being the active site, i.e. where 
sites n and n + 1 are represented without truncation (cf. 
Eq. 0): This way, the basis states chosen to represent 
optimally the state before Un is applied, 

= i^aar0M'j)n\T)n+l\/3) (7) 

are equally well suited for representing the state 

Unl'lp) = ^{Un)aT,a'T'i^aa'T'/3\(^)\cr)n\r)n+l\P) (8) 

ac7T f3 

a' t' 

without any additional error, because Un only acts on 
the part of the Hilbert space (spanned by the vectors 

n\T) n+i) that is exactly represented. 

To continue the sweep, a DMRG truncation is car- 
ried out with Un\'^) being the target state instead of |?/^). 
The key observation is that the new truncated basis is 
optimally adapted to t/n|V^) and different from the one 
that would have been chosen for |?/^). In contrast to the 
conventional static DMRG^^, the optimally represented 
Hilbert space hence follows the time- evolution of the state 

Then basis transformations to the left or right are per- 
formed, until the next part of Eq. can be applied. We 
thus apply the full operator of Eq. © by sweeping the 
active site n through the system. The price to be paid 
is that a truncation error is introduced at each iteration 
step of the sweep as is known from static DMRG. 

To start time-dependent DMRG, some initial state has 
to be prepared. There is no unique recipe, the most ef- 
fective one depending on the desired initial state. The 
procedure we adopt for our initial state |ini) is to calcu- 
late it as the ground state of a suitably chosen Hamilto- 
nian i^ini (which does in principle not have to have any 
physical significance). Such a choice is i^ini = BnS^^ 
with Bn < for n on the left, Bn > for n on the 
right half of the chain. In this case, a physical picture for 
i^ini does exist; it corresponds to switching on a magnetic 
field that aligns the spins and that is strong enough for 
all interactions in Eq. Q to be negligible. 



IV. ACCURACY OF THE ADAPTIVE 
TIME-DEPENDENT DMRG 

As so far no quantitative analysis of the accuracy of 
the adaptive t-DMRG has been given in the literature. 



we provide a detailed error analysis for the time evolution 
of the initial state |ini) in a spin- 1/2 quantum XX chain, 
i.e. Jz = S = 0. This system is an excellent benchmark 
for the adaptive t-DMRG due to its exact solution^ that 
can be compared to the DMRG results. The exact so- 
lution reveals a nontrivial behaviour with a complicated 
substructure in the magnetization profile. From a DMRG 
point of view this Hamiltonian is not too specific in the 
sense that the experience from static DMRG suggests a 
relatively weak truncation error dependence on J^. 



A. Possible errors 

Two main sources of error occur in the adaptive t- 
DMRG: 

(i) The Trotter error due to the Trotter decomposi- 
tion. For an nth-order Trotter decomposition^^, the error 
made in one time step dt is of order dt^~^^. To reach a 
given time t one has to perform t/dt time-steps, such 
that in the worst case the error grows linearly in time t 
and the resulting error is of order {dt)^t. In our setup 
of the Trotter decomposition, the error scales linearly 
with system size L, and overall it is of order {dt)^Lt for 
the times of interest. (Eventually, the error must satu- 
rate at a finite value, as measured quantities are typically 
bounded.) The linear L dependence of the error is ex- 
pected for generic initial states. For the particular choice 
of |ini) of this paper, however, many of the 0{L) contri- 
butions to the Trotter error vanish, as many of the sites 
exhibit no dynamics at all for short times. - For the cal- 
culations presented in this paper, we have chosen n = 2, 
but our observations should be generic. 

(ii) The DMRG truncation error due to the representa- 
tion of the time-evolving quantum state in reduced (al- 
beit "optimally" chosen) Hilbert spaces and to the re- 
peated transformations between different truncated basis 
sets. While the truncation error e that sets the scale of 
the error of the wave function and operators is typically 
very small, here it will strongly accumulate as 0{Lt/dt) 
truncations are carried out up to time t. This is because 
the truncated DMRG wave function has norm less than 
one and is renormalized at each truncation by a factor 
of (1 — e)~^ > 1. Truncation errors should therefore 
accumulate roughly exponentially with an exponent of 
eLt/dt^ such that eventually the adaptive t-DMRG will 
break down at too long times. The error measure we 
use here saturates at 0(1) and sets a limit on the expo- 
nential growth; also, partial compensations of errors in 
observables may slow down the error growth. The ac- 
cumulated truncation error should decrease considerably 
with an increasing number of kept DMRG states m. For 
a fixed time t, it should decrease as the Trotter time step 
dt is increased, as the number of truncations decreases 
with the number of time steps t/dt. 

At this point, it is worthwhile to mention that our 
subsequent error analysis should also be pertinent to the 
very closely related time-evolution algorithm introduced 
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by Verstraete et a/.^^, which differs from ours for the 
present purpose in one major point: In our algorithm a 
basis truncation is performed after each local apphcation 
of Un- In their algorithm truncations are performed af- 
ter all local time-evolutions have been carried out, i.e. 
after a global time-evolution using U = fin 
iterative procedure, the wave function after such a full 
time evolution is not guaranteed to be the globally op- 
timal state representing the time-evolved state. How- 
ever, for small dt the state update via the operators Un 
is likely to be small, we expect the global optimum to be 
rather well approximated using the present algorithm, as 
seems to be borne out by direct comparisons between 
both approaches'^. Errors should therefore exhibit very 
similar behaviour. 

We remind that no error is encountered in the applica- 
tion of the local time evolution operator Un to the state 
lip)^ as is discussed after Eq. ((Hj). 



B. Error analysis for the XX-model 

In this section, we analyze the errors from the adaptive 
t-DMRG in the time evolution of the XX-model by com- 
paring it to the exact solution^, with the ultimate goal of 
finding optimal DMRG control parameters to minimize 
the errors. 

We use two main measures for the error: 

(i) As a measure for the overall error we consider the mag- 
netization deviation the maximum deviation of the local 
magnetization found by DMRG from the exact result, 

err(t) = max^l (^^,dmrg(^)) - (^n,exact(^))|. (9) 

In the present study, the maximum was typically found 
close to the center of the chain. 

(ii) As a measure which excludes the Trotter error we 
use the forth-back deviation FB{t)^ which we define as 
the deviation between the initial state |ini) and the state 
\fb{t)) = U {—t)U {t)\mi) ^ i.e. the state obtained by evolv- 
ing |ini) to some time t and then back to t = again. If 
we Trotter-decompose the time evolution operator U{—t) 
into odd and even bonds in the reverse order of the de- 
composition of U{t)^ the identity U{—t) = U{t)~^ holds 
without any Trotter error, and the forth-back deviation 
has the appealing property to capture the truncation er- 
ror only. In contrast to the magnetization deviation, the 
forth-back error does not rely on the existence of an exact 
solution. 

As our DMRG setup does not allow easy access to the 
fidelity | (ini|/6(t)) |, we define the forth-back deviation to 
be the L2 measure for the difference of the magnetization 
profiles of |ini) and \fb{t)), 

FBit)=(Y,i(^m\S^Jmi)-{fb{t)\S^Jfbit))A . (10) 
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FIG. 3: Magnetization deviation as a function of time for 
different Trotter time steps dt and for m = 50 DMRG states. 
At small times (region A in the inset), the deviation is dom- 
inated by the linearly growing Trotter error for small times. 
At later times (region B in the inset), much faster, non- linear 
growth of the deviation sets in at some well-defined runaway- 
time tR. As shown in the inset, tn increases with increasing 
dt. 
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FIG. 4: Magnetization deviation as a function of Trotter time 
step dt (system size L — 100, m = 50 DMRG states) at times 
t = 5 (upper figure) and t = 30 (lower figure). For t = 5, the 
magnetization deviation is quadratic in dt as expected from 
the Trotter error. For t — 30, at small dt the magnetization 
deviation is no longer quadratic in dt and larger than the Trot- 
ter error would suggests. This is a signal of the contribution 
of the accumulated truncation error. 
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FIG. 5: The forth-back error FB{t) for t = 30 and t = 50, 
as function of dt. Here, L — 100, m — 50. 



In order to control Trotter and truncation error, two 
DMRG control parameters are available, the number of 
DMRG states m and the Trotter time step dt. 

To study the effect of varying dt^ consider the mag- 
netization deviation as shown in Fig. |S1 Two main ob- 
servations can be made. At small times (regime A), the 
magnetization deviation decreases with dt and is linear 
in t as expected from the Trotter error. Indeed, as shown 
in the upper part of Fig. ^ the magnetization deviation 
depends quadratically on dt for fixed t, and the Trotter 
error dominates over the truncation error. At large times 
(regime B), the magnetization deviation is no longer lin- 
ear in t, but grows almost exponentially, and also does no 
longer show simple monotonic behaviour in dt: The mag- 
netization deviation in this regime is obviously no longer 
dominated by the Trotter error, but by the accumulated 
truncation error. 

The two regimes A and B are very clearly separated 
by some runaway time tR^ with regime A for t < tR and 
regime B for t > tR (a precise procedure for its deter- 
mination will be outlined below). The runaway time tR 
increases when dt is increased: Because the total num- 
ber of Trotter time steps t/dt is decreased, the accumu- 
lated truncation error decreases, and the Trotter error in- 
creases, hence the competing two errors break even later. 

This (it-dependence of tR is also seen in the lower part 
of Fig. ^ where the dt dependence of the magnetization 
deviation is plotted at some larger time {t = 30) than 
in the upper part, t = 30 is larger than the runaway 
time (i.e. in regime B) for dt < 0.05, in regime A other- 
wise. We see indeed for dt > 0.05 (region A) the familiar 
quadratic Trotter error dependence. For small dt < 0.05 
(region B), the deviation is dominated by the accumu- 
lated truncation error that increases as dt decreases. This 
is reflected in the growth of the magnetization deviation 
as dt is decreased. 

The almost exponential growth of the truncation error 
with the number of Trotter steps can also be seen from 
the forth-back deviation that is not susceptible to the 
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FIG. 6: Magnetization deviation AM(t) as a function of 
time for different numbers m of DMRG states. The Trotter 
time interval is fixed at dt = 0.05. Again, two regimes can 
be distinguished: For early times, for which the Trotter error 
dominates, the error is slowly growing (essentially linearly) 
and independent of m (regime A); for later times, the error is 
entirely given by the truncation error, which is m-dependent 
and growing fast (almost exponential up to some saturation; 
regime B). The transition between the two regimes occurs at 
a well-defined "runaway time" tn (small squares). The inset 
shows a monotonic, roughly linear dependence of tn on m. 



Trotter error. In Fig. O we show the forth-back devia- 
tion FB{t) for t = 30 and t = 50 as a function of the 
Trotter time step dt. FB(t) increases as a consequence 
of the stronger accumulation of the truncation error with 
decreasing Trotter step size dt and hence an increasing 
number of steps t/dt. 

Let us now consider the dependence of the magnetiza- 
tion deviation err(t) on the second control parameter, the 
number m of DMRG states. In Fig.lHl err(t) is plotted for 
a fixed Trotter time step dt = 0.05 and different values 
of m. In agreement with our previous observations, some 
m-dependent "runaway time" tR^ separates two regimes: 
for t < tR (regime A) , the deviation grows essentially lin- 
early in time and is independent of m, for t > tR (regime 
B), it suddenly starts to grow more rapidly than any 
power-law. The onset of a significant m-dependence has 
indeed been our operational definition of tR in Fig. Eland 
ini In the inset of Fig. El tR is seen to increase roughly 
linearly with growing m. As m ^ oo corresponds to 
the complete absence of the truncation error, the m- 
independent bottom curve of Fig. El is a measure for the 
deviation due to the Trotter error alone and the runaway 
time can be read off very precisely as the moment in time 
when the truncation error starts to dominate. 

That the crossover from a dominating Trotter error at 
short times and a dominating truncation error at long 
times is so sharp may seem surprising at first, but can be 
explained easily by observing that the Trotter error grows 
only linearly in time, but the accumulated truncation er- 
ror grows almost exponentially in time. The latter fact is 
shown in Fig. [3 where the forth-back deviation FB(t) is 
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FIG. 7: The forth-back error FB{t) for L = 100, m = 40, 
dt — 0.01 and dt — 0.05, as function of t. 
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FIG. 8: The forth-back error FB(t) for t = 50 and t = 30, 
as function of m. Here, L = 100, dt = 0.05. 



plotted as a function of t for some fixed m. Here, we find 
that the effects of the truncation error are below machine 
precision for t < 10 and then grow almost exponentially 
in time up to some saturation. 

By comparison, consider Fig. [HI where FB{t) is plot- 
ted as a function of m, for t = 30 and t = 50. An 
approximately exponential increase of the accuracy of 
the method with growing m is observed for a fixed time. 
Our numerical results that indicate a roughly linear time- 
dependence of tR on m (inset of Fig. ^ are the conse- 
quence of some balancing of very fast growth of precision 
with m and decay of precision with t. 

Before concluding this section, let us briefly consider 
a number of other possible effects that might affect 
tR. One might alternatively conceive that the well- 
defined runaway-time tR results from a sudden failure 
(of stochastic or of fundamental nature) of the trunca- 
tion algorithm to capture one important basis state. It 
can be refuted on the basis of Fig. [51 Fig. [7j and Fig. [SI 
Such an error should manifest itself as a pronounced step 
in FB(t)^ depending on the time evolution having gone 
past tR or not. Such a step is, however, not observed. 



XX model: entropy of half chain 
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time 



FIG. 9: Entanglement entropy Se from Eq. (|11|) between the 
left and the right half of the chain as function of time. 

tR might also be thought to reflect a fundamental 
DMRG limit, namely a growth of the entanglement 
within the time-evolved state which the limited number 
of DMRG states m is not able to capture adequately at 
t > tR. This scenario can be excluded by observing the 
strong dependence of tR on the number of time steps, 
which this scenario cannot explain. Indeed, a study of 
the entanglement entropy between the left and the right 
half of the chain 

Se(t) = Trplog^p, (11) 

p being the reduced density matrix of the left (or equiva- 
lently the right) half of the chain, confirms this view: As 
shown in Fig. [31 Se{t) is only mildly growing with time 
and is well below the maximum entanglement entropy 
5'max ~ log2^ ^hat the DMRG can reproduce. 

Therefore we conclude that the error at short times 
is dominated by the Trotter error, which is independent 
of m and approximately growing linearly with time. At 
some runaway time, we observe a sharp crossover to a 
regime in which the m-dependent and almost expontially 
growing truncation error is dominating. This crossover is 
sharp due to drastically different growth of the two types 
of errors. The runaway time thus indicates an immi- 
nent breakdown of the method and is a good, albeit very 
conservative measure of available simulation times. We 
expect the above error analysis for the adaptive t-DMRG 
to be generic for other models. The truncation error will 
remain also in approaches that dispose of the Trotter er- 
ror; maximally reachable simulation times should there- 
fore be roughly the same or somewhat shorter if other 
approximations enhance the truncation error. 

C. Optimal choice of DMRG parameters 

How can the overall error - which we found to be a del- 
icate balance between the Trotter and the accumulated 
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truncation error - be minimized and the important run- 
away time be found in practice? From the above scenario 
it should be expected that the truncated density matrix 
weight at each step does not behave differently before or 
after the runaway time and hence is no immediately use- 
ful indicator to identify the runaway time. This can in 
fact be seen from Fig. where the truncated weight is 
shown for the same parameters as in Fig. IHl Also, it is 
not obvious to extract a precise relationship between the 
truncation errors at each DMRG truncation and the ac- 
cumulated errors. Instead, a precise convergence analysis 
in m or dt seems to be more telling and easily feasible. 

Of course, it is desirable to choose the number of kept 
states m as large as possible within the constraints re- 
garding the available computer resources. This choice 
having been made, the runaway time tR is determined 
for different Trotter time steps dt by comparing different 
values of m as in Fig. El Only two slightly different val- 
ues of m are sufficient for that purpose. Now the Trotter 
time step dt is chosen such that the desired time t is just 
below tR. This way, the optimal balance between the 
Trotter error and the truncation error is found, which 
corresponds in the lower part of Fig. ^ to the minimum 
of err(t) on the border between regime A and B: The to- 
tal error would increase at larger dt due to the Trotter 
error, and at smaller dt due to the truncation error. 

Thus, it is a good practice to choose for small times 
rather small values of dt in order to minimize the Trotter 
error; for large times, it makes sense to choose a some- 
what coarser time interval, in order to push the runaway 
time to as large values as possible. 

In terms of numbers of time steps, we conclude from 
Fig. [31 that for the present model and our parameters 
(L = 100 - 200), the adaptive time-dependent DMRG 
seems to be able to perform about 1000-5000 time steps 
reliably even for m = 50, depending on the desired level 
of accuracy, corresponding to 0(100/ J) in "real" time. 
We note that this is a very small value of m by DMRG 
standards, and that using an optimized code, one should 
be able to increase m by an order of magnitude, and 
hence access much longer times (by an order of magni- 
tude). 



V. LONG-TIME PROPERTIES OF THE 
TIME-EVOLUTION 

the time evolution of the initial state |ini) on 
the XX chain at temperature T = was examined in the 
long-time limit using the exact solution. It was found 
that the magnetization Sz{n^t) given in Eq. Q can be 
described for long times in terms of a simple scaling func- 
tion, Sz{n^t) ~ ^ {{n — nc)/t)^ where ric is the position 
of the chain center. The scaling function is the solution of 
the partial differential equation dtSz + dxj{Sz) = with 
the magnetization current j{Sz) = I/ttcos IttS'^I which 
has been shown to describe the macroscopic time evolu- 
tion of the magnetization profilei. The characteristics. 
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FIG. 10: The lost weight in the density matrix trunca- 
tion, summed over time intervals At = 0.1, is shown for the 
same parameters as in Fig. |3 A comparison with Fig. |3 re- 
veals, however, that both values are not useful criteria for the 
DMRG truncation error and are in particular not suited to 
reveal the runaway time tR. 

i.e. the lines of constant magnetization S'^, have a slope 
V = sin Itt^^I- 

The magnetization profile ^ {{n — nc)/t) has a well- 
defined front at (n — nc)/t = ±1, i.e. is moving outwards 
ballistically with velocity v = 1. On top of this overall 
scaling form an additional step-like substructure arises, 
which was analysed in detail in 21]. It was found that 
while the step width broadens as t^/^, the step height 
decreases as t~^/^, such that the integrated transported 
magnetization within each step remains constant at 1. It 
was suggested that each of these steps corresponds to a 
localized flipped spin flowing outwards. 

The XX model, however, has several very special prop- 
erties: It corresponds to a free-fermion model and is 
therefore exactly solvable; it is critical; and its total cur- 
rent operator J = ^^jn commutes with the Hamilto- 
nian, [J, i7] = 0. One may ask whether the above find- 
ings are due to any of the particularities of the XX model 
or more generic. 

The adaptive t-DMRG allows us to study the long-time 
evolution of |ini) in different coupling regimes of Eq. 
We chose two extensions of the XX model, namely a 
S^S^- interaction, and dimerization. 

In Fig. ^2 and El we visualize the time evolution of 
the local magnetization in density plots, with site index 
n on the x-axis, time t on the ?/-axis. Here, the absolute 
value of the magnetization is shown as a grayscale and in 
lines of constant magnetization at [(S'^)! = 0.2, 0.4. In 
Fig. ^2 the relation between the density plots and the 
actual magnetization profile for the XX model is shown 
at two times, t = and t = 40. The exact solution is 
perfectly reproduced, including the detailed substructure 
of the magnetization profile. 

In Fig. El density plots for various values of Jz be- 
tween and 1.1 are shown. For small Jz {Jz < 1), we 
observe ballistic transport of the magnetization. This 
regime is characterized by a constant transport velocity 
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FIG. 11: Left: Time evolution of the absolute value of the 
local magnetization \{Sn{t))\ for the XX model as a density 
plot, where the local magnetization itself is exactly antisym- 
metric with regard to the chain center. The lines of constant- 
magnetization (Sn) = ±0.2, dbO.4 are shown as solid lines. As 
an illustration, local magnetizations (5'^(t)) for the time slices 
t = and t = 40 are shown explicitly. A step-like substruc- 
ture can be seen for t = 40 in perfect quantitative agreement 
with the exact solution. Error bars are below visibility. 



of the magnetization, hence the lines of constant magne- 
tization shown in Fig. E| are approximately straight for 
Jz < 1. The magnetization front propagation slows down 
as Jz increases, and almost comes to a halt when > 1. 
Although the sharpness of this crossover dX Jz = 1 is sur- 
prising, its general nature can be understood from the 
limits Jz ^ ^ and \Jz\ oo: For small ^ the 
S^S^- and 5'^S'^-interactions dominate. Being spin flip 
terms, they smear out the initially hard step profile in 
the z magnetization. For large J^, on the other hand, 
the 5'^5'^-interaction dominates. This term does not de- 
localize the step profile, and in the limit \Jz\ oo, the 
initial state is even a stationary eigenstate of the Hamil- 
tonian. 

Besides the structure of the overall front, we also ob- 
serve for Jz ^ ^ remnants of the steplike substructure 
from the XX model, individual pockets of transported 
magnetization at velocity 1, which we call "carriers". As 
Jz is increased, these carriers keep the velocity v ~ 1, but 
are increasingly damped and thus less and less effective 
in transporting magnetization. 

In order to put the above observations on a more quan- 
titative footing, we plot in Fig. El the integrated flow of 
magnetization through the center. 



AM(t) = / {jL/2{t'))dt' 
Jo 



L 

E 

2>L/2 



{{S^{t)) + l/2). (12) 



This quantity has the advantage that unlike the lines of 
constant magnetization in Fig. and El it shows the 
overall spin transport without being too much biased by 
single "carriers". 

We observe in Fig. El roughly linear behaviour of 
AM{t) for \Jz\ < 1, which suggests ballistic magneti- 
zation transport at least on the time scales accessible to 
us. As Jz increases, magnetization transport slows down 



until around Jz = I the behaviour changes drastically: 
For > 1, AM{t) seems to saturate at a finite value, 
around which it oscillates. On the time scales accessi- 
ble to us, we thus find a sharp crossover at = 1 from 
ballistic transport to an almost constant magnetization. 

This crossover is even more clearly visible in Fig. 
where we plot the exponent a of the magnetization, 
AM{t) (X t^, for values Jz between and 1.5. Here, 
the exponent a is close to 1 for < 1, confirming the 
roughly linear transport, and quickly drops to zero in the 
regime of constant magnetization for > 1. 

Fig.Elillustrates how the exponent a was obtained, for 
the special case = 1. Here the exponent a = 0.6 ± 0.1 
indicates that the magnetization transport is clearly not 
ballistic anymore. In fact, we find from a scaling plot 
Fig. El that for long times the magnetization collapses 
best for a scaling function of the form Sz{n^ t) ~ (j){n/t^-^) 
with an uncertainty in the exponent of approximately 
0.1, indicating superdiffusive or diffusive transport in the 
time range under consideration. 

The proposed crossover from ballistic to almost no 
transport is also visible in the expectation value of the 
current jn = Jnl^{{S:^ S~^iO)). For Jz = S = 0, it is 
known^ that the current at the middle of the chain ap- 
proaches a finite value as t ^ oo. This is only possible 
for ballistic transport. In the case of (sub- / super-) dif- 
fusive transport or constant / oscillatory magnetization, 
on the other hand, the central current must fall off to 
zero as the magnetization gradient flattens or must even 
become negative to allow for the oscillations. 

This expected behaviour is seen in Fig. El where we 
plot the current at the center of the chain as a function 
of time for various values of Jz between and 1.1. We 
averaged the current over the 5 middle sites in order to 
filter out local current oscillations. We observe that for 
relatively long times, the current approaches a constant 
value for \Jz\ < 1, whereas the current falls off rapidly 
and then seems to exhibit damped oscillations around 
zero for \Jz\ > 1. This strengthens our previous conclu- 
sion of a crossover from ballistic transport to a more or 
less constant magnetization at \Jz\ = 1. 

Remarkably, this crossover for the behaviour of a high- 
energy quantum state |ini) is found at the location Jz = I 
of the quantum phase transition from the critical phase 
to the Neel antiferromagnetic state (see Fig. a priori 
a low-energy event. To understand the subtle connection 
between the time evolution of |ini) and the phase transi- 
tion, we exploit that the time-evolution does not depend 
on the sign of Jz^ as discussed in Sec. ^ Therefore the 
time evolution of the high-energy state |ini) for Jz > I 
is identical to that for = — Jz < —1, where |ini) is a 
low-energy state. At the quantum phase transition from 
the ferromagnetic state to the critical phase at = — 1 
the ground state, a kink state for < — 1 (if we impose 
the boundary condition spin up on the left boundary and 
spin down on the right boundary) changes drastically 
to a state with no kink and power-law correlations for 
> —1. Therefore, our initial state is very close to an 
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FIG. 12: Density plots of the magnetization \{Sn{t))\ as in Fig. El the values of Jz being (from left to right, top to bottom) 
0, 0.3, 0.6, 0.9, 1.0, 1.1, and ^ = 0. For better visibility of the profile, the grayscale mapping of \{Sn{t)) \ was chosen differently 
in each plot as indicated by the legends. Solid lines: lines of constant magnetization (S^) = ±0.2, ±0.4; these allow for a direct 
comparison of the magnetization transport between different Jz- The ray- like structure indicates the "carriers". 
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FIG. 13: The change in the magnetization AM(t) is 
shown. The curves are plotted in the order Jz = 
0; 0.3; 0.6; 0.9; 1.0; 1.1; 1.5, where = is the steep- 
est. The curves Jz = 0; 0.3; 0.6; 0.9 show the same linear 
behaviour for the observed times, i.e. up to t = 60. 



eigenstate - the ground state - for < —1, but not for 
> —1. Thus, the harsh change in the time-evolution 

of the high-energy state |ini) at = 1 can be explained 

by the severe change in the ground state properties at 
= — 1 , and the crossover is linked to a quantum phase 

transition at a different location in the phase diagram. 
We now study the influence of a nonzero dimerization 



FIG. 14: Best fit for the exponent a in AM(t) oc t"", for 
the data shown in Fig. 1131 and for times between t = 20 and 
t = 60. We estimate the uncertainty in a to be of the order 
of 0.1 due to the limited time available (cf. Fig. I15f) . It was 
not possible to fit the slow oscillations for = 1.1. To the 
eye, however, the curve in Fig. suggests slow oscillations 
around a constant value, hence we included in the data point 
a = for J2 = 1.1 by hand (encircled). 



S in Eq. We restrict our analysis to the case Jz = 0. 
The dimerized models can still be described in terms 
of the free-fermion picture and are exactly solvable (for 
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FIG. 15: = 1: The change of the magnetization in a 

double logarithmic plot with an algebraic fit. 
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FIG. 16: = 1: Collapse of magnetization for a superdiffu- 
sive scaling form (x/t*^"^). 




FIG. 17: Current, averaged over the 5 middle sites, for vari- 
ous values of Jz between and 1.1. 



static properties see 0). The current, however, is not 
conserved for nonzero dimerization. This example will 
shed light on the question whether the long-time limit 
depends on current conservation or on the free-fermion 
property, or yet on other special properties of the system. 
We expect two obvious effects of nonzero dimerization: 
Firstly, the overall front velocity should slow down, be- 
cause the magnetization now propagates faster on half of 
the links, but slower on the other half, the net effect be- 
ing a reduction of the total velocity. Secondly, we expect 
oscillations with a period of two lattice sites. This is ob- 
vious in the limit ^ ^ 1, where each strongly coupled pair 
of sites can be viewed as an almost isolated subsystem, 
in which the magnetization oscillates back and forth. We 
expect remnants of this behaviour also at dimerizations 
\5\<1. 

The data shown in Fig. confirms this expectation 
qualitatively, but does not reveal any other qualitative 
change of the long-time limit for nonzero dimerization. 
For = 1, the system is trivially given by isolated pairs 
of neighbouring sites, therefore the propagation velocity 
drops to zero. 

Fig. ^1 and Fig. [201 reveals explicitly that no qualita- 
tive change occurs as the dimerization is switched on: 
the change in magnetization AM{t) still shows the lin- 
ear behaviour typical of ballistic transport. For increas- 
ing S ^ 1 oscillations on top of this linear behaviour 
arise. We find that switching on finite dimerization does 
not change the long-time behaviour of the time evolution 
also for nonzero Jz (not shown). In particular, the time 
evolution here is drastically influenced by the transition 
at = 1 as in the case ^ = discussed above. 

To summarize, we find the same long-time behaviour 
of the initial state | j ... jj ... |) in the dimerized sys- 
tem — a system with gapped excitation spectrum and 
which is exactly solvable — as in the system with small 
^^^^-interaction, \Jz\ < 1 — a system which is critical 
— whereas the behaviour changes drastically for larger 
6'^6'^-interaction, \Jz\ > 1. Hence we cannot attribute 
the ballistic transport of the magnetization to the spe- 
cific properties of the XX model; neither to be exactly 
solvable, nor to the continuous spectrum nor to the con- 
served current in the XX-model. The drastic change at 
\Jz\ = 1 stems from the special property of the initial 
state to resemble the ground state in the ferromagnetic 
phase and the highest energy state in the antiferromag- 
netic phase. 

Finally, let us include a note on the errors in the 
present analysis. A convergence analysis in m as in sec- 
tion 1^ shows that the errors and the runaway time are 
roughly the same as for the XX model. The plot in 
Fig. El goes up to time t = 95, whereas the runaway 
time tR is somewhat earlier, ^ 60 — 80, depending 
on the precise value of Jz- Indeed, a convergence anal- 
ysis in m reveals that the accuracy in the central region 
decreases for t > Ir. For dimerized models the run- 
away time tR is somewhat shorter (between Ir = 40 and 
tR = 80 for m = 50, depending on the dimerization). 
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FIG. 18: Density plots of the magnetization {Sn{t)) as in Fig. for dimerization (from left to right, top to bottom) 
(5 = 0; 0.2; 0.4; 0.6; 0.8; 1.0, and Jz = 0. The grayscale mapping is different in each plot as indicated by the legends. Solid lines: 
lines of constant magnetization (Sn) — ±0.2, ±0.4. 
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FIG. 19: Change in magnetization AM(t) for different dimer- 
izations, from top to bottom: ^ = 0, 0.2, 0.4, 0.6, 0.8, 1.0. 



FIG. 20: Best fit for the exponent a in AM(t) oc t"", for 
the data shown in Fig. 1131 and for times between t = 20 and 
t = 40. 



This fact reflects the reduced accuracy of the DMRG al- 
gorithm when dealing with inhomogenous systems. As 
always, it is possible to increase tR by increasing m. 



VI. CONCLUSIONS 

We have investigated the evolution of the initial state 
I t ... 11 ... I) under the effect of nearest-neighbour 



interactions with the adaptive time-dependent DMRG. 

For weak 6'^6'^-interaction, i.e. IJ^I < 1 in Eq. (QJ, 
and arbitrary dimerization, < S < 1, we find that for 
long times the transport of the magnetization is ballistic 
as it was found for the XX-model^. The magnetization 
profile shows the same scaling form for long times, i.e. 
S'^(n, t) = (p{{n — nc)/t), where ric is the position of the 
chain center, but with different scaling functions cp. For 
stronger 6'^5'^-interaction, i.e. \Jz\ > 1, even in a ho- 
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mogeneous system, ^ = 0, a drastic change in the long- 
time evolution is seen. The magnetization transport is no 
longer ballistic, but shows oscillatory behaviour around 
a constant value. Hence our results suggest that the spe- 
cific properties of the XX model are not responsible for 
ballistic transport at long times. The drastic change in 
the long time behaviour at the phase transition = 1 
can be attributed to the close resemblance of the initial 
state to the ground state for < — 1. 

Our error analysis for the adaptive time-dependent 
DMRG shows that for small times the error is dominated 
by the Trotter error whereas for long times the truncation 
error becomes the most important. This finding should 



be general and hold for non-exactly solvable models as 
well, and should therefore allow to control the accuracy 
of the results of adaptive time- dependent DMRG in gen- 
eral models. Overall, we find this method to be very 
precise at relatively long times. 
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